The A7V — > NN weak interaction in effective field theory 
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The nonleptonic weak \JAS\ = 1 AN interaction, responsible for the dominant nonmesonic decay 
of all but the lightest hypernuclei, is studied in the framework of an effective field theory. The long- 
range physics is described through tree- level exchange of the SU(3) Goldstone bosons while the short- 
range potential is parametrized in terms of lowest-order contact terms. We obtain reasonable fits to 
available weak hypernuclear decay rates and quote the values for the parity-violating asymmetry as 
predicted by the present effective field theory. 



PACS numbers: 21.80.+a, 13.30.Eg, 13.75.Ev 



For the past 50 years, A- hypernuclei, systems of one 
or more A hyperons bound to a core nucleus, have been 
used to extend our knowledge of both the strong and the 
weak baryon-baryon interaction from the NN case into 
the SU(3) sector. The nonmesonic hypernuclear weak de- 
cay, facilitated via the \AS\ — 1 four-fermion interaction, 
thus complements the weak AS* = NN case, which al- 
lows the study of only the parity-violating amplitudes. 

In analogy to JVAf phenomenology, the nonmesonic hy- 
pernuclear decay has traditionally been modeled using a 
meson-exchange approach [1] . The long-range part of the 
interaction is naturally explained by one-pion exchange 
which could approximately reproduce the total (one- 
nucleon induced) nonmesonic decay rate, AN — > NN, 
but not the partial rates, the proton-induced Ap — > np 
rate T p and the neutron- induced An — > nn rate r n [2]. 
Due to the A7V mass difference, the process AN — > NN 
produces nucleons with momenta around « 420 MeV, 
suggesting that the short-range part of the interaction 
cannot be neglected. These contributions have been de- 
scribed either through the exchange of vector mesons 
[3,4], whose production thresholds are too high for the 
free A decay, or through direct quark exchange [5]. 

In contrast to previous theoretical studies, we present 
an exploratory study in order to determine the possi- 
ble efficacy of the use of Effective Field Theory (EFT) 
methods in hypernuclear decay. Studies in this direction 
have already begun in Ref. [6] where a Fermi (V-A) in- 
teraction was added to the OPE mechanism to describe 
the weak AN — > JVJV transition. The present approach 
is motivated by the remarkable success of EFT tech- 
niques based on chiral expansions in the (non-strange) 
SU(2) sector [7-10], which suggests its extension to the 
SU(3) realm, even though stability of the chiral expan- 
sion is less clear for the SU(3) sector, due to the signifi- 
cant degree of SU(3) symmetry breaking. A well-known 
example of the problems facing SU(3) chiral perturba- 
tion theory has been the prediction [11-14] of the four 



parity-conserving (PC) amplitudes in the weak nonlep- 
tonic decays of octet baryons, Y — > Nir, with Y = A, S 
or 2. In particular, Refs. [14] and [15] studied the con- 
tributions from negative-parity intermediate states and 
demonstrated their potential to resolve this longstand- 
ing issue. For our purpose, these higher-order effects are 
beyond the simple lowest-order analysis considered here. 

The EFT approach is based on the existence of well 
separated scales in the physical process under study. For- 
mally, the high-momentum (short-distance) modes in the 
four-baryon interaction Lagrangian are replaced by con- 
tact operators, compatible with the underlying physical 
symmetries, of increasing dimension. Built in such way, 
the Lagrangian will contain an infinite number of terms, 
and a consistent power-counting scheme is needed in or- 
der to truncate such an expansion to a given order. The 
coefficients appearing in the Lagrangian are then fitted 
to reproduce the available data in the low-energy regime. 
Whether an EFT will succeed to describe a particular 
process or not is directly related to the success in ob- 
taining a controlled, systematic expansion in terms of a 
small parameter. In contrast to the NN case, however, 
the AN — > JVJV transition corresponds to an energy re- 
lease w 177 MeV (|p| m 417 MeV) at threshold. It is 
therefore not at all clear if low-energy expansions can 
be successfully carried out. In light of this energy re- 
lease at threshold, it is reasonable to include the pion 
(m* sa 138 MeV) and the kaon (m K « 494 MeV) as 
dynamical fields. Working within SU(3) also supports 
treating the pion and kaon on equal footing. The last 
member of the SU(3) Goldstone-boson octet, the 77, is 
usually not included, since the strong rjNN coupling is 
an order of magnitude smaller than the strong ttNN and 
KAN couplings [16]. Thus, following a power-counting 
scheme based on the engineering dimensions of the oper- 
ators, at leading order (order unity in the external mo- 
mentum, p° ) the present study of the weak AN — > NN 
transition includes the contribution of the long-ranged 
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pion and kaon exchanges, and a short-range contribution 
given by leading nonderivative contact terms (leading or- 
der parity-conserving, LO PC, terms). Our study also 
includes next-to-leading order (NLO) terms in the mo- 
mentum expansion (or equivalcntly, leading order parity- 
violating pieces, LO PV). Not considered here are con- 
tributions from intermediate-range 27r-exchange. Such a 
component is two orders higher in the chiral expansion 
than the corresponding single pion-exchange piece. Also, 
previous studies [17] found such contributions to be small 
due to significant cancellations between the correlated 
and uncorrelated 27r-exchanges. 

The weak and strong Lagrangians for pion and kaon 
exchanges are the same as in Ref. [18], and we use exper- 
imental values for the couplings at the strong (g NN7T = 
13.16) and weak (ANir = 1.05 for the parity-violating 
amplitude and —7.15 for the parity-conserving one, in 
units of Gpfn^ 2 = 2.21 x 1CP 7 ) baryon-baryon-pion ver- 
tices. Like the NNtt coupling, the ANK and YjNK cou- 
pling constants represent a fundamental input into our 
calculation. Unlike the NNir coupling, however, their 
values are considerably less well known, with <7ank — 
13 — 17 and g SNK = 3 — 6. Here, we choose the val- 
ues given by the Nijmegen Soft Core NSC97f interaction 
model [19],g ANK = —17.66 and g SNK = —5.38. Of course, 
the weak nuclcon-nucleon-kaon coupling constants are 
not accessible experimentally, so we obtain numerical 
values by making use of SU(3) and chiral algebra con- 
siderations [3,4]. In addition, those potentials will be 
regularized by using monopole form factors at each ver- 
tex [20]. We note that all the strong model-dependent 
ingredients used in the present calculation (such as cut- 
off parameters or strong couplings) have been taken from 
the Nijmegen Soft-Core model f [19]. 

If no model is assumed, the low energy A7V — ► NN 
process can be parametrized through the 6 partial waves 
listed in Table I [22]. The % -» % and 3 Si -> 3 5i 
transitions can only be produced by the 1 • S 3 (f) and 
(?i • c?2 • <5 3 (r ) operators, where 5 3 (f ) represents the con- 
tact interaction. The 1 S'o — ► 3 Po and 3 Si — > x Pi tran- 
sitions proceed through the combination of the spin- 
nonconserving operators: (a\ — 02) ■ {pi — p>2 , 8 3 {f) }, 
(<7i -<r 2 ) • [pi -p 2 , S 3 (r)}, i(<?i xct 2 )-{pi-P2, 8 3 (r)} 
and i (t?i x B2) • [pi — P2 , S 3 (f) ], where pi is the derivative 
operator acting on the "itft" particle [23]. The 3 Si — > 3 P\ 
transition is allowed by the combination of the spin- 
conserving operators [B\ + 0*2) • {pi — Vi , 8 3 {f) } and 
(c?i + 1? 2 ) • [pi — P2 , 5 3 (r ) ], while only two-derivative op- 
erators can produce the last (tensor) transition. 

Using power counting we can discard operators of order 
q 2 /Mat 2 , where the momentum transferred is defined by 
q = pa — PN t — Pzv 2 — Pn- The remaining (lowest order) 
operators lead to the following four-fcrmion interaction 
in r-space (in units of G F = 1.166 x lO" 11 McV" 2 ): 
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where the last factor represents the AI = 1/2 isospin 
part of the 4P interaction. Note that the delta functions 
have been smeared by using a Gaussian form with a typi- 
cal vector-meson (p) range, 5 <~ \Plm~ 1 w 0.36 fm. Here, 
M = (3M + Ma)/4 is a weighted average of TV, A masses 
while C s and C 3 P are the ]th Low Energy Coefficients 
(LEC) at 0th and first order respectively. Although the 
form of the contact terms is model- independent, the size 
of these LEC's depends upon how the theory is formu- 
lated, and they are expected to be of order of the other 
couplings in the problem. These couplings provide a very 
simple representation of the short distance contributions 
to the process at hand. In a complete model, they would 
be represented by specific dynamical contributions, such 
as p, lj, etc. -exchange. However, we eschew the tempta- 
tion to be more specific — in fact this generality is one 
of the strengths of our approach. We evaluate the coef- 
ficients purely phenomenologically and leave theoretical 
interpretation of the pieces to future investigations. Of 
course, the specific size of such coefficients depends upon 
the chiral order to which we are working. However, if 
the expansion is convergent, then the values of these ef- 
fective couplings should be relatively stable as NLO or 
higher effects are included. 

It is well known that the high momentum transferred 
in the AiV — ► NN reaction makes this process sensi- 
tive to the short range physics which is characterized by 
our contact coefficients. Moreover, since the |A5| = 1 
reaction takes place in a finite nucleus, the extraction 
of reliable information of the elementary weak two-body 
interaction requires a careful investigation of the many- 
body nuclear effects present in the hypernucleus. In the 
present calculation, we use a shell-model for the initial 
hypernucleus, where the single-particle A and N orbits 
are taken to be solutions of harmonic oscillator mean 
field potentials with parameters (6a and &n) adjusted to 
experimental separation energies and charge form factor 
(respectively) of the hypernucleus under study. For A 2 C 
and 6 N = 1.64 fm and b A = 1.87 fm, while for ^He 
they take the values &n = 1-4 fm and 6a = 1-85 fm. The 
strong YN interaction at short distances, absent in mean- 
field models, is accounted for by replacing the mean-field 
two-particle AA^ wave function by a correlated [28] one 
obtained from a microscopic finite-nucleus G-matrix cal- 
culation [29] using the soft-core and hard-core Nijmegen 
models [30] . The NN wave function is obtained from the 
Lippmann-Schwinger (T-matrix) equation with the input 
of the Nijmegen Soft Core potential model; details of the 
calculation can be found in the Appendix of Ref. [18] 
which presented a detailed study of different approaches 
to Final State Interactions (FSI) in the decay process. 
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We begin the discussion of our results with a remark 
on the data. One might wonder if there can be only 
three independent data points in the nonmesonic decay: 
the proton-induced and neutron-induced rates T p and 
r„, and the asymmetry A (associated with the proton- 
induced decay), relating observables from one hypernu- 
cleus to another through hypcrnuclcar structure coeffi- 
cients. While one may indeed expect measurements from 
different p-shell hypernuclei, say, A=12 and 16, to pro- 
vide the same constraint, the situation is different when 
including data from s-shcll hypernuclei like A=5. For 
the latter, the initial AN pair can only be in a relative 
s-state, while for the former, relative p-states are allowed 
as well. We therefore include data from the A=5,ll and 
12 hypernuclei in our fits. 

Note that we do not attempt to perform a final, quan- 
titative fit to all hypcrnuclear data. Rather, we are ex- 
ploring whether an EFT can be used to reproduce various 
reasonable subsets of the hypcrnuclear decay data in or- 
der to verify the validity of such an expansion. We note 
that the presently large experimental error bars in hy- 
pernuclear decay observables puts strong limitations in 
any EFT approach to the decay process. In this sense, 
the values for the parameters we are presenting have to 
be taken with caution. A quantitatively more rigorous 
understanding awaits better data as well as theory which 
includes higher-order effects. Only the more recent mea- 
surements from the last 12 years were used, excluding, 
however, those r„/r p (hereafter = n/p) ratio data whose 
error bars were larger than 100%. A more detailed dis- 
cussion of the minimization calculation will be presented 
elsewhere [24,25]. 

We also have taken the data at face value and have not 
applied any corrections due to, e.g., the two-nucleon in- 
duced mechanism, which has been estimated to amount 
up to 25% of the total decay rate for p-shell hypernu- 
clei [26,27]. Ideally, exclusive experiments would separate 
this mechanism from the measured total decay rate, per- 
mitting a fit to observables that are not contaminated 
with multi-nucleon effects. At the present time, given 
the sizable error bars of the data, this omission would 
not change our conclusions. 

No parameters were fitted for the results with only 
7r and K exchange, shown in Table III. As has been 
known for a long time, tt exchange alone reasonably 
well describes the observed total rates, while dramat- 
ically underestimating the n/p ratio. The tensor PC 
channel dominates the proton-induced rate while it is ab- 
sent in the L = neutron-induced one. Incorporation of 



kaon exchange yields a destructive interference between 
both mechanisms (OPE and OKE) in the PC amplitudes, 
while the interference is constructive in their PV coun- 
terparts. As a consequence, n/p is enhanced by about 
a factor of five, within reach of the lower bounds of the 
experimental measurements, while the total rate under- 
predicts the observed value by about a factor of two. This 
-k — K interference also leads to values for the asymmetry 
that are close to experiment for the p-shell hypernuclei, 
but far off for A=5 . Since the contributions of both 
ry-exchange and two-pion exchange are negligible, these 
discrepancies illustrate the need for short-range physics. 

Allowing contact terms of order unity (leading-order 
PC operators) to contribute leads to four free parameters, 
Cg, C s , Cis and Civ- Data on the total and partial de- 
cay rates for all three hypcrnuclear systems are included 
in the fit, but no asymmetry measurements. The inclu- 
sion of the contact terms roughly doubles the values for 
the total decay rates, thus restoring agreement with ex- 
periment. The impact on the n/p ratio is noteworthy: 
the value for ^He increases by 10% while the n/p ratios 
for and \ 2 C almost double. This is an example of 
the differing impact certain operators can have for s- and 
p-shell hypernuclei. The effect on the asymmetry is op- 
posite, almost no change for A=ll and 12, but a 30% 
change for A=5. The magnitudes of the parameters, Cg, 
Cg, Civ, listed in Table II, are each around their natural 
size of unity, while Cis is a factor of three or so larger. 
Note the substantial error bars on all the parameters, 
reflecting the uncertainties in the measurements. 

Three new parameters are admitted when we allow the 
leading-order PV terms (of order q/M^) to contribute 
with the coefficients C P , Cp, and Cp. As shown in Ta- 
ble II, the parameters for the PV contact terms are larger 
than the ones for the PC terms, and in fact, compatible 
with zero. Including the three new parameters does not 
substantially alter the previously fitted ones, thus sup- 
porting the validity of our expansion. Regarding their 
impact on the observables, the PV contact terms barely 
modify the total and partial rates but significantly affect 
the asymmetry, as one would expect for an observable de- 
fined by the interference between PV and PC amplitudes. 
The calculated asymmetry changes sign for all three hy- 
pernuclei, moving the ^He value within the measured 
range at the expense of the one for "B. This shift oc- 
curs without any asymmetry data constraining the fit. In 
order to further understand this behavior, we have per- 
formed a number of fits including the asymmetry data 
points of either ^He or "B or both. Tables II and III 



1 Note that for A He we quote the value of the intrinsic A asymmetry parameter, a^, which is experimentally accessible, while for 
p-shell hypernuclei the accessible quantity is A, the difference between the number of protons coming parallel and antiparallel 
with respect to the polarization axis. This quantity can be related to a a through the relation A = paoa, where pa is the A 
polarization, to be extracted from theoretical models. 
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display the result of one of those fits. Inclusion of the 
^He asymmetry helps in constraining the values of two 
of the LO PV parameters. We find that the two present 
experimental values for A=5 and A=ll cannot be fitted 
simultaneously with this set of contact terms. Future 
experiments will have to settle this issue. 

We have also performed fits allowing a contribution 
from an isospin AI = 3/2 transition operator. The 
resulting fit, with \ 2 ~ 1.4, shifts strength from the 
isoscalar contribution to the new A J = 3/2 one, leaving 
the other parameters unchanged. However, as shown in 
Table II, we can clearly get an good fit to all observables 
without such transitions while obtaining couplings of rea- 
sonable size. In addition, we have checked that our con- 
clusions are independent of the strong interaction model 
used to describe FSI in the transition. Employing NN 
wave functions that are obtained with either the NSC97f 
or the NSC97a model in the fit leaves the observables 
almost unchanged, with the exception of the asymmetry 
parameter, which can change up to 50%. The obtained 
couplings can easily absorb the changes but remain com- 
patible within their error bars. Similarly, we performed 
a study of the sensitivity of the calculated observables to 
the smearing function in Eq. (1). For values of <5 going 
from 0.3 to 0.4 fm, the results are remarkably insensitive, 
except again in the case of the asymmetry. 

In conclusion, we have studied the nonmesonic weak 
decay using an Effective Field Theory framework for the 
weak interaction. The long-range components were de- 
scribed with pion and kaon exchange, while the short- 
range part is parametrized in leading-order PV and PC 
contact terms. We find coefficients of natural size with 
significant error bars, reflecting the level of experimental 
uncertainty. The largest contact term corresponds to an 
isoscalar, spin-independent central operator. There is no 
indication of any contact terms violating the AI =1/2 
rule. In this study we have not speculated as to the dy- 
namical origin of these contact contributions. Rather our 



aim was to ascertain their basic magnitude and to estab- 
lish the validity of the EFT framework for the weak decay. 
The next generation of data from recent high-precision 
weak decay experiments currently under analysis holds 
the promise to provide much improved constraints for 
studies of this nature. 
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TABLE II. LEC coefficients corresponding to the LO cal- 
culation. The values in parentheses have been obtained in- 
cluding aA(iH e ) in the fit. 





+ LO (PC) 


+ LO (PC±PV) 




-1.54 ±0.39 


-1.31 ±0.41 (-1.04 ±0.33) 


ch 


-0.87 ±0.24 


-0.70 ±0.35 (-0.57 ±0.27) 


C° P 




-5.82 ±5.31 (-4.49 ±1.57) 






2.47 ±3.13 (1.84 ±1.93) 


cl 




-5.68 ±3.13 (-4.47 ±2.31) 


Cis 


5.01 ± 1.26 


4.68 ±0.67 (5.97 ±0.86) 


Civ 


1.45 ±0.38 


1.22 ±0.20 (1.56 ±0.26) 



TABLE III. Results obtained for the weak decay observables, when a fit to the T and n/p for \Be, A *B and is performed. 
The values in parentheses have been obtained including a^Qjie) in the fit. 
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± LO 
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0.42 
0.09 
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0.23 
0.50 
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0.43 
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0.55 (0.55) 
0.28 (0.24) 
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1.04tg;™ [31] 
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0.08 
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0.41 
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-0.06 


0.95 
0.67 
-0.05 


0.93 (0.93) 
0.77 (0.77) 
0.03 (0.02) 


1.14 ±0.2 [31], 0.89 ±0.15 [32], 0.83 ± 0.11 [36] 
0.87 ± 0.23 [38] 
-0.01 ±0.10 [35] 
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0.93 


1.54 (1.15) 
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